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TECHNICAL NOTE NO. 1636 

DETERMINATION OF STRESSES IN GAS-TURBINE DISKS 
SUBJECTED TO PLASTIC FLOW AND CREEP 
By M. B. Millenson and S. S. Mans on 


SUMMARI 

A finite-difference method previously presented for computing 
elastic stresses in rotating disks is extended to include the com- 
putation of the disk stresses when plastic f lov and creep are con- 
sidered. A finite-difference method is employed to eliminate numer- 
ical integration and to permit nontechnical personnel to make the 
calculations with a minimum of engineering supervision. Illustra- 
tive examples are included to facilitate explanation of the pro- 
cedure by carrying out the computations on a typical gas -turbine 
disk through a complete running cycle. 

The results of the numerical examples presented indicate 
qualitatively that plastic flow markedly alters the elastic -stress 
distribution and that, if the amount of creep Is small, the effect 
on stress distribution is also small. 


INTRODUCTION 

With the advent of jet propulsion as a motive force for air- 
craft, the gas turbine has become an important source of power. In 
most machinery, design stresses are limited by the yield strength 
or the creep strength of the material employed, together with a 
certain factor of safety, and little or no analytical consideration 
is given to the occurrence of plastic flow under operating condi- 
tions. Gas-turbine disks, however, are required to operate under 
thermal gradients and centrifugal forces producing stresses that, 
in materials currently available, frequently exceed the yield 
strength, resulting in plastic flow. The interaction cf plastic 
flow and creep, together with the variation of thermal gradients 
through a series of cycles consisting in starting, running, and 
stopping can produce stress distributions and even failures that 
might not be suspected on a basis of elastic -stress analysis. 

A rapid routine method cf elastic -stress analysis of rotating 
disks is presented in reference 1, which gives accurate values of 
the true stresses in disks, provided that the yield strength of the 
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material la not exceeded. The finite-difference method, cf refer- 
ence 1 has been extended at the HACA Cleveland laboratory to 
include consideration of - plastic flow and creep, which thus allows 
calculation of the true stresses in a gas-turbine disk and gives 
the variation of stress distribution with time. The handling cf 
plastic flow is somewhat less routine than the calculation of the 
elastic stresses in that a repetitive trial procedure is required. 
With practice, the correct value can be obtained on the fourth or 
fifth trial. The computation of the effect cf creep, although in 
procedure the same as the computation cf plastic flow, is a direct 
calculation requiring no tidal -and -error procedures. Because the 
method eliminates numerical integration, nontechnical personnel can 
make the calculations with a minimum of engineering supervision. 


SYMBOLS 

The following symbols are used: 
c creep rate under stress o Q , inches per inch per hour 
E elastic modulus of disk material, pounds per square inch 
h axial thickness cf disk, inches 

B ratio 


r radial distance, inches 
T temperature, °P 
u radial displacement, inches 

a coefficient of thermal expansion between act ual temperature 
and temperature at which there is zero thermal stress. 
Inches per inch per °F 

T total creep under stress ct q , inches per inch 

A plastic increment of strain, inches per inch 

& r plastic Increment of strain in radial direction 

Ajj plastic increment of strain in tangential direction 

AT temperature increment above temperature of zero thermal 
stress, 
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8 r creep increment in radial direction, inches per inch 
8-j. creep increment in tangential direction, inches per inch 

€ strain, inches per inch 

€ plastic strain corresponding to stress a e in tensile spec! 
men, inches per inch 

e r radial strain, inches per inch 

tangential strain, inches per inch 

p Poisson’s ratio 

O A 

p mass density of disk material, pound second per inch* 

a stress, pounds per square inch 

o Q equivalent tensile stress, pounds per square inch 

O r radial stress, pounds per square inch 

Oj. tangential stress, pounds per square inch 

a proportional elastic limit, pounds per square inch 

y 

T time during which creep occurs, hours » 

co angular velocity, radians -per second 

The following supplementary subscripts are used for denoting 
values of the preceding symbols in connection with the finite- 
difference solution: 

n n^* 1 point station 
"bll 

n-1 (n-1) point station 

a station at smallest disk radius considered (Far disk with a 
central hole, this station is taken at the radius of hole; 
for a solid disk this station is taken at a radius approxi 
mately 5 percent of the rim radius.) 

b station at rim of disk or base of blades 
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The following supplementary symbols denote combinations of the 
foregoing symbols: 


Ar,n 

■H,n 

®r,n 

®t,n 


> 




stress coefficients defined by equations 
a r,n = -^r^n a t,a + ®r,n 
a t,n = A t,n °t,a + ®t,n 


C n = r n *n 


^n (! + » n )(* n - r^) 
» " ^ + 2 ^ r n 

^n = | ( r n - r n-l) *n 


_1_ d + Mn)(r n - r a . a ) 
n ** + 2 ®n r n 

F n c r n-l Vl 

y, Vl q ♦ ^n-l )(r n - r n-l> 

n " Vl 2 E n-1 r n-l 

= 2 <*n “ r n-l^ ^n-1 

_1_ . q + ^n-lH r n - r n-l) 
u * E n-1 2 Vl r n-l 

*n = f w2 ( r n - r n-l)(Pn ^n r n 2 + Pn-1 Vl r n-l ) 



*n 


<*n AT n " “n-l AT n-l 

F n ^n ~ F n D n 
C'n^a -ChD’h 


v 
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K» 


C n F» n - C' n F n 
n - C* a T± - C n D« a 


a'n^ + Gh D 'n 
c ’n D n " c n D ’n 


L ’n' 5 


^ - 


M' 


C*n ^ + C n S» n 
" c n D n - c n D 'n 

V'nZn + Vn fr'n ~ p n ~ Q'n) 

C* D - C D* 
n. n n n 

C»n *n + On ( H 'n - p> n ~ <*'n) 
c 'n D n - c n D 'n 


(M^ and M' n are defined in reference 1 for the special oase 


P'n - Q' a = 0) 


% ” r n-l\ 


( T n ’ r n-l\ ( 

J r,n \ 2 r n ) + A r,n-1 y 2 r n-1 j " A + ~ |1 + 


u,n 


r n ~ r n- l 
2 r_ 


+ A t,n-1 1 - T7 


r n ‘ r n-l> 


n-1 


Q' 


8 , 


Tr, - r, 


n-1 


r,n i 2 r_ 


"W Ft^t) - \ n (l - ^ 


r_ - r. 


n-1' 


2 r„ 


r n - r n-^ 


+ ^t,n-l 1 " 2 r, 


n-1 


MALTS IS CF PLASTIC FLOW AHD CHEEP 

Assumptions . - Four assumptions are made in the subsequent 
analysis : 

1. The dish material is linearly elastic up to a limiting stress 
value, called the proportional elastic limit, and above this limit 
plastic flow occurs. 
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2. All variables cf material properties and operating condi- 
tions are symmetrical about the axis of rotation. 

3. Axial stresses may be neglected and the radial and tan- 
gential stresses are uniform across the thickness cf the disk. 

4. Temperatures are uniform, across the thickness cf the disk. 

Outline of method. - In any thin, rotating disk, the complete 
stress state is defined •when the two principal stresses, radial a r 
and tangential a t , are known at every radius . Two equations 

relating these stresses to the radius are required to specify the 
stress distribution. The first of these equations can be deter- 
mined from the conditions cf equilibrium of an element of the disk 
and involves no elastic properties of the material. The second is 
derived from the compatibility conditions, which state the inter- 
relation of radial and tangential strains. The compatibility 
conditions are dependent upon stress -strain phenomena and must 
therefore include any departure from linear elasticity. When 
modification to allow for any possible departure from Hooke’s law 
is made, the compatibility conditions become t rue for any value of 
stress. The equation derived from the compatibility conditions 
thus modified, together with the equilibrium equation, is treated 
by the finite -difference method of reference 1, and similar equa- 
tions are obtained. These equations result in additional terms in 
the final equations, which are used to modify the result of the 
elastic calculation. 

Whenever stresses under discussion have been calculated by the 
method of reference 1 only, they will be referred to as "elastic 
stresses" j where plastic flow and creep have been taken into 
account, the stresses will be referred to as "plastic stresses." 

Derivation of method. - The equilibrium equation, which applies 
to both the elastic and plastic conditions, is 

~ (r h a r ) - h a t + p a> 2 r 2 h = 0 (l) 

The elastic compatibility relations given in terms of the radial 
displacement are 


du a r ~ ° t 


E 


+ a AT 


( 2 ) 
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and 


6fc 


u 

r 



+ a AT 


( 3 ) 


Equations (2) and (3) must "be modified to include consideration of 
plastic flow. When a material is stressed "beyond the proportional 
elastic limit, the strain in the material is different from that 
indicated "by Hooke's law. The strain under such a load may "be 
considered as "being made up of two components, one elastic as pre- 
dicted "by the laws of elasticity and one an increment of strain due 
to the flow that oocurs. Bewriting equations (2) and (3) on this 
"basis gives 

H . % + a M + ^ (5) 


Similarly, any creep that occurs represents an additional 
departure from elastic behavior, which further modifies equa- 
tions (2) and (3) to 


du °r ~ » a t 

€ r = 3? = J + a AT + Ay + 5 r 


u 

•t “ F = 


°+. - H 


E 


+ a AT + A^ + &tj 


( 6 ) 

( 7 ) 


When the parameter u is eliminated as in reference 1, 


1 

E Q t 


-|a P + aAT + A+ 




1-Hl 


A _A . 




Applying the finite -difference method to equations (l) and (8) 
and using the notation introduced in the section entitled "SXMBOIS” 
gives 


°n a r,n " ^ °t,n “ F n a r,n-l + Gp a t,n-l -B n 


( 9 ) 
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and 


' a „ 
n r,n 


- D 


'n a t,n 


F 'n a r,n-l " ^'n °t,n-l + H 'n “ B 'n " 


(10) 


The solution of the equations is facilitated hy the substitu- 
tion of the stress coefficients A^j v A^ n , ®r,ni and- ®t,n 

into equations (9) and (10). Proceeding as in reference 1 results 
in the equations 


On ^,n“-®*n ^r>Q-l”^i ^t>n-l - 0 

O’n A^ n -F’ n ^ n _i^' n ^n-l - 0 

c n B r,n"^u B t,n“ F n ^^n-l'On B^n-l+Bn = 0 
0'n B r,n'’ I), n B t,n -F 'n ^r^n-l+^’n ^t^n-l^'n+^'n+^’n = 0 



All hut the last of equations (ll) and equations (15) of refer- 
ence 1 are identical. When equations (ll) are solved for Aj. n , 
•H^n* ®r,n* ani Bt,n> 


A r,n = ^n-l + In A t,n-1 
A t,n = K 'n Ar,n-1 + L ’n ^tjn-l 
®r,n = K n ^r^n-l + In ^t^n-l + 

®t,n = K 'n ®r,u-l + l*n Bt,n-1 + M* n 


( 12 ) 


The symbols K^, K' n , In* aa( ^ ^’n kave same meaning 
as in reference 1. The and M' n ' terms are now defined as 




p, n ^ + Pn ( B> n - ?'n 
C ’n D n - C n D 'n 



( 13 ) 
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M’ n = 


c 'n*n + c n fr'n “ p ’n ~ Q'n> 


n* tl - n_ -n'. 


(13a) 


The elastic case of reference 1 thus "becomes a special case of the 
more general problem in which P' n and Q' n are both zero. 


Evaluation of plastic terms. - In order to apply the finite- 
difference method to problems involving plastic flow, a relation 
between stresses and strains in the plastic region must be esta- 
blished. In reference 2, a numerical-integration method for 
computing dish stresses is presented in which elongation is assumed 
to proceed at constant stress when the proportional limit is reached. 
References 3 and 4 present equations for the plastic relation of 
stress to strain based on the maximum distortion theory. Relations 
can be derived from the equations given in reference 4, which form 
a convenient means of finding the plastic increments corresponding 
to the stresses present in the disk. Rewriting these equations in 
the notation of this report and letting the subscripts 1, 2, and 3 
denote the three principal directions in the most general case gives 


A 1 = f [(<*. " ° 2 ) + (<*i - <7 3 >] 

d 2 = f [(02 - Ox) + (a 2 - a 3 )] 

A 3 = | [° 3 - a l> + < a 3 - <>2>] 


a e = ^(oi - o 2 ) 2 + (a 2 - a 3 ) 2 + (a 3 - a x ) 2 


(14) 


(15) 


where the ratio R is defined in terms of the corresponding uniaxial 
stress <J Q and plastic strain Cp in a tensile speoimen by the 
relation 


R 



(16) 


When equations (14), (15), and (16) are reduced to the biaxial con- 
dition, which is assumed to prevail in the disk (that is, a 3 * 0), 
and the finite-difference notation is introduced 
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A r,n - ^ ( 9 r,n " a t,n + a r,n) 

A t,n “= 3 ( a t,n " °r,n + a t,n^ 

a e,n = ^°r, n 2 “ a r,n °t,n + a t,n 2 


1 


and 


3 £ 

E = &£ 


2 a 


e.n 


Substituting equation (19) into equations (17) gives 


k r,n “ 2 <j®|“ °r,n - °t,n) 
4 t,n * (2 ' a r,n) 


'e.n 


(17) 


(18) 


(19) 


( 20 ) 


A typical uniaxial stress-strain curve illustrating the relation 
between effective stress a 0jI1 and effective plastic strain £p,n 
on such a curve is shown in figure 1. Investigations at the National 
Physical Laboratory of Great Britain on turbine-disk alloys and 
experiments by Taylor and Quinney (reference 5) were found to 
correlate well with equations ( 20 ) for the range of strains over 
which the volume of the material is approximately constant. 


Equations (20) give the relations between plastic strains and 
true stresses that will be used as the basis for numerical calcula- 
tions in the present report. The method of stress analysis to be 
presented does not depend, however, on the validity of these equa- 
tions. As more accurate relations are determined between stresses 
and strains, these relations may readily be used in place of 
equations ( 20 ) . 


Calculation of plastic flow when no previous plastic flow has 
occurred. - The determination of the plastic stresses in the disk 
resolves itself into the problem of finding corresponding stresses 
and strains that satisfy equilibrium and compatibility equations (9) 
and (10), and biaxial stress-strain equations (20) . The problem is 
approached by first computing the elastic stresses, and the equiva- 
lent uniaxial tensile stress at each station is determined from 
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equation (18). If at any station this stress exceeds the propor- 
tional elastic limit of the material at the temperature at this 
station, then plastic flow takes place, and it "becomes necessary to 
resort to a trial -and-error procedure to adjust the stresses to 
allow for this flow. 

Assume, for example, the equivalent uniaxial stress at a given 
station lies at point A on the extension of the modulus line in 
figure 2. Because the point A lies above the proportional elastic 
limit (point B), plastic flow must occur. The stress and strain 
must "be adjusted to fall on the curved stress-strain curve that is 
characteristic of the material. As a starting point, the total 
strain in the true stress-strain condition is assumed equal to the 
strain at A. The stress-strain condition at the given station then 
lies on the constant-strain line through B, or at C. The plastio 
strain £p^ n is given by CD. Values of Aj^ n , A t ^ n , and P' n 
may be obtained by using this value of £ p ^ n , together with the 
values of (J r ^ n , cr t,n^ an ^ °e,n ** oni the elastic calculations. 

Once P* n has been calculated, new values of a r>n , < 7 -j^ n , 
and a e ^ n can be computed. The new value of Q e ,n ie greater 
than that at point D, such as that at point E. Although the stresses 
corresponding to G e ,n at point E together with the strain CD meet 
the conditions of equations (9) and (10), they locate the stress- 
strain point F, which is not on the stress-strain curve, so that 
the physical conditions Imposed by the material are as yet unsatis- 
fied. Inasmuch as any value of €p^ n lass than CD would give a 
value of cr e ^ n greater than that at E, CD Is a lower limit 
of Cp n . Similarly, because the value of a e ^ n calculated by 
using an €p^ n of CD is too great, the increment of strain EG 
corresponding to this 0" e ^ n is an upper limit of €p,n* Inasmuch 
as the true value of £p >n lies between CD and EG, their numerical 
average, shown as HK, is assumed to be a good approximation. Hew 
values of P’ n , <?r,n, a t,n> a-ud a e,n 03X1 computed by using 
HK for Cp^, the stress at E for C e ^ n , and a r ^ n and 
Suppose this new value of O e ^ n lies at the point M. Because the 
stress at M is higher than the stress at H in value, the increment HK. 
Is too small a value for €p^ n and is therefore established as a 
new lower limit of £p,n* Further, because M is less than E, the 
corresponding increment M is a new upper limit for €p jtl and the 
process could be repeated again with the numerical average of MN 
and HK. Similarly, if the calculation using an €p jn of HK had 
resulted in a <J e ,n a-fc F, HK would constitute a new upper limit 
and PQ a new lower limit. Had the resulting a Qj n ^ een at B, HK would 
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have still became the new upper limit of £p,n> CD would have 

remained as the lower limit. The process is repeated until the value 
of Cp n used in the computation and the Cp^ n corresponding to the 
resulting a e ^ n are equal. 

Calculation of plastic flow when previous plastic flow has 
occurred . - The equations for strain that would apply to a disk that 
had already undergone the plastic strain are 

£ r = + oAT + [Ay] + Ay (21) 

35 

€ t = + oAT + [A t 3 + At (21a) 

Here the terms [Ay] and [At] represent strains already 
existent in the material before the application of stresses at 
and Gy and are constant for the calculation, whereas A r and At 
represent the components of plastic strain resulting from the 
application of G r and Gf In the solution of the equations by 
the finite -difference method, a term [P' n ] appears together with 
term P' n . When previous plastic. flow has occurred only once, 

Cp ' n 3 is identical with P’ n from the previous calculation; where 
plastic flow has previously occurred more than once, [P’ n ] is the 
algebraic sum of all earlier P' n terms. ThuB, the previous plastic 
flow given by [P* n ] may be grouped with the temperature-effect 
tern H' n by replacing H' n with H' n - [P' n ]. 

This procedure amounts to an assumption that, as the load and 
the temperature change, the stress position on the new stress-strain 
curve would be the same as if a test specimen were loaded above the 
yield point, the load removed, the temperature changed, and a new 
load applied. This assumption is illustrated by figure 3, in which 
point A represents a loading at. the first temperature condition; the 
dotted line AB represents the load-removal path; the curve BCD, the 
stress-strain curve at the new temperature; and point C, the new 
stress position. The total strain at this point C is given by the 
sum of three strains. The residual strain caused by the first 
loading is e-jj eg is the elastic part of the strain caused by the 
second loading; and € 3 , the plastic strain caused by the second 
loading. 

When the foregoing procedure is applied, the curve BCD must, 
of oourse, represent the true stress-strain curve at the new tem- 
perature of a material that has already been subjected to the plastic 
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cycle OAB. In general, the new stress-strain curve is different 
from the stress-strain curve at the given temperature of a material 
that has not been subjected to plastic flow; however, unless data 
are available it may be necessary to assume that the curve BCD is 
the stress-strain curve at the given temperature of a spec imen of 
virgin material. 


Calculation of effect of creep. - Creep is usually defined as 
the continuous deformation of material under a continuously applied 
load. Experimental data on creep of various materials are usually 
obtained from tests run under constant load and temperature, aitfrnHgh 
in many engineering applications of materials the more general 
problem of changing load and temperature must be considered. The 
deformation curve obtained in a typical test is shown in figure 4. 
From this figure it can be seen that the deformation may be 
considered as having occurred in three stages. During the primary 
stage, the deformation proceeds at a decreasing rate; during the 
secondary stage, at a constant rate; and during the tertiary stage, 
at an increasing rate, which proceeds until failure occurs. 


Because of the lack of data on creep except for nnlA-rinl 
tensile stress, a relation between creep deformation and stress 
must be assumed. The following equations have been used for cal- 
culations in this report but, as better data become available, 
more accurate relations can be used. By the use of reasoning 
similar to that employed in determining the biaxial components of 
plastic-strain formulas for the creep increments. By n and 8-^ n 
may be written * ’ 


8 


r,n 


n 


2 a. 


< 2 °r,n 



( 22 ) 



2 a, 


e,n 



(22a) 


In equations (22) and (22a), F n represents the total creep 
that would occur in time T under the uniaxial stress o e n . It 
Is here assumed that for sufficiently small values of T £he creep 
may be considered as ooourring instantaneously at the end of the 
time period. 


During the secondary stage of creep, a characteristic creep 
rate Cq exists, corresponding to the stress <j e n a t temper- 
ature T, and T n is given directly by 


r n ° °n T 


(23) 
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This rate is the -value usually published in papers on creep and 
is the rate used for the numerical calculations of this report. 
During primary and tertiary creep stages, the creep rate is alBo a 
function of time, hut does not otherwise complicate the computation. 

Once -values of 5 r ^ n and 5t,n have been found, the values of 
the Q' n terms may be determined and new values of Or, n ®nd Ot,n 
may be computed. If the computed values of a r ^ n and CT.^ n . differ 
by more than a small amount, perhaps 2 percent, from the values of 
these stresses before creep occurred, a shorter time interval should 
be selected and additional computations made for eaoh such time 
interval required to equal the total time during which creep occurs. 
The effect of creep that occurred at previous time intervals is 
considered in a manner similar to that employed in considering 
previous plastic flow. The successive values of Q' n are summed to 
form a term [Q' n 3, which gives the total effect of all previous 
creep deformation so that the term 

H'n - CP‘ n 3 


is replaced by 


H'n - tP’n] - 

In any calculation of stress distribution subsequent to the 
occurrence of creep, the creep term [Q' n ] is combined with the 
term [P' n ] as the cumulative effect of all previous plastic 
deformation. 

Examples showing in detail how successive stages of plastic flow 
and creep are computed, each stage considering all previous plastic 
deformation, are given in the section entitled "NUMERICAL EXAMPLES." 


NUMERICAL EXAMPLES 

The numerical examples presented here represent a set of com- 
putations during one complete start -run- at op cycle for a typical 
turbine disk with a continuous rim and welded blades. The assumed 
profile of the disk is shown in figure 5, together with the loca- 
tions of the point stations used in the computations. The assumed 
temperature distributions and corresponding turbine rotative 
speeds are shown in figure 6. Curve IT and the corresponding speed 
of 11,500 rpm represent the steady-state -r unning condition. 

Curves I to III and the corresponding speeds represent running 
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conditions through which the turbine disk passes in reaching steady- 
state operation. Curves Y to Vli together with the respective 
speeds represent running conditions through which the turbine disk 
passes when being stopped. Creep is assumed to occur only during 
the steady-state running period. 

The physical properties of the disk material, including specific 
gravity, modulus of elasticity, stress-strain characteristics, and 
thermal coefficients of expansion, were based on the data appearing 
in references 6 and 7, together with unpublished data obtained from 
the author of these references. The stress-strain curves were 
constructed on the basis of these data and those for example I 
appear in figure 7. Inasmuch as no data were available on the 
effect of previous plastic flow on the shape of the stress-strain 
curves, it was necessary to ignore such effects and to use curves 
obtained directly from simple tensile-test data. The creep 
properties used in the computations were assumed, but they corres- 
pond approximately to the secondary creep rates given in reference 8. 
The effect of primary creep was emitted because of lack of data. 

Because the disk used for these calculations is solid at the 
center, a supplementary numerical example showing the Computation 
of the plastic-flow effect on stress distribution in a disk 
containing a central hole is given in the appendix. 

Example I, - Example I is the calculation of the stress dis- 
tribution in a disk operating under the conditions of curve I of 
figure 6 and having been subjected to no previous plastic deforma- 
tion. These conditions sire assumed to represent disk operation 
after the first short period of steady combustion, when gas tem- 
peratures are high, thereby establishing a steep temperature 
gradient between the center and the rim of the disk. 

The preliminary elastic calculation is carried out in 
table I (a) by the method of reference 1. Two changes are made in 
the tabular setup. The first change is the insertion of 
columns 25a and 25b immediately following column 25. Column 25a * 
lists the accumulated values of [P' n ] and • [Q’ n ], the total 
effect of previous plastic deformation. For the present example, 
this col umn is zero for all stations. Column 25b is the value of 
the term H' n - [P' n ] - [Q' n ], which in this example is the same 
as column 25. The second change is the computation of Ma and M' n 
(columns 31 and 32), which were computed in reference 1 by the use 
of column 25. In the present computations, column 25b is used. In 
addition, two more columns, 40 and 41, are added. Column 40 lists 
the values of the proportional elastic limit Oj f n of the material 
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and column 41 lists the values of G e n as computed from equa- 
tion (18). The entries in columns 40* and 41 of table l(a) show that 
the equivalent stress C e ^ n is less than Oy^ n for all point 
stations except 17 to h. The effect of plastic flow must he 
considered at these stations and flow at these stations modifies the 
stresses at other locations in the disk. With the exceptions and the 
additions noted, the method of computation is the same as the method 
of reference 1 and will not he discussed in further detail. 

The plastic-flow calculation has been divided into two parts 
because several quantities used in the computation depend only on 
the dimensions of the disk and can be used in all subsequent calcu- 
lations involving plastic deformation. These quantities are 
computed for stations 17 to b as shown by the four column headings 
of table 1(b). 

The second part of the plastic-flow calculation is given in 
table 1(c) . The first column in this part of the table (column 46) 
lists the values of e p ^ n obtained from the corresponding stress- 
strain curve (fig. 7) as explained previously. Column 46a lists 
the value of e p ^ n used for the ensuing calculation, which for the 
first approximation is the same as column 46. Columns 47 and 43 
list the values of A r ^ n and A^ n , respectively, computed by 
equations (20) . Columns 49 to 52 are computed as shown by the 
column headings and from these columns the values of P' n are 
computed and listed in column 53. Col umn 54 gives the values of 
the term H' n - [P T n ] - [Q'J - P'n, which is then used to compute 
new values of M’ n , B^ n , Bt >n , G^ n , O., 

and q^ n as shown in columns 55 to 62, respectively. The new 
values of € Pjn corresponding to the new values of a e ^ n are read 
from figure 7 and listed in column 46 of the second-approximation 
calculation. The values in column 46 for the first and second 
approximations now constitute the lower and upper limits, respectively, 
of the possible strain increments. For the second approximation, 
column 46a therefore lists as the values of e p ^ n to be used in this 
set of calculations the numerical averages of the two sets of 
readings from the stress-strain curve. Fran this- value, another new 
set of stres® values is computed and a third set of readings listed 
in column 46. 

At this point in the calculation, two alternate procedures are 
possible, as shown by consideration of station b. Inasmuch as the 
average value of 4300 x 10 -6 inches per inch used in the second 
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approximation gave a graph reading of 1900 x 10 inches per inch, 
the averaging procedure would indicate that the next trial should he 

- t . - 5 . 9 - 6 .- x 10" 6 or 4130 x 10 -6 inches per inch 

This value could be used and the procedure continued until the correct 
value is found. Considerable time may be saved in making the calcu- 
lation, however, ‘if a weighted approximation is used. Because the 
plastio- strain value of 3960 x 10“ 6 inches per inch gave a resulting 
reading of 4650 x 10"® inches per inch whereas the value 4300 x 10“° 
inches per inch gave the reading 1900 x 10"® inches per inch, the 
strain 3960 x 10”® inches per inch is* apparently more nearly correct 
than 4300 x 10"® inches per inch. In addition, the shape of the 
stress-strain curve in the region of 3960 x 10"® is suoh that small 
increases in stress correspond to large changes in strain. If a 
trial calculation were made using a value closer to 3960 then 4130 
(for instance, 4100), more information might be obtained than would 
be obtained by the averaging procedure. The right answer is thereby 
obtained more quickly. The sauna reasoning might be applied to the 
selection of values to be used at the other stations for the third 
calculation. The second of these two procedures has been used in 
table 1(c), as can be seen from the values of €p jn In column 46a 
used for the third approximation. 

Completion of the third approximation and comparison with new 
values of strain obtained from the stress-strain curve shows the 
estimates of the third approximation to be nearly correct, so that 
small adjustments made to compute the fourth and fifth approximations 
give the final answers. A calculation equi valent to a sixth approxi- 
mation is then made to column 53 to get the final correct values of 
the P ' n terms. The stresses at the other stations a through 16 
can now be computed by using the value of <Jt,a found in the sixth 
approximation with the values A r ^ n , A-t^n, Br >n , and n foun d 
in table 1(a). The values of plastic stress at all radii together 
with the elastic-stress distribution are plotted in figure 8. 

Example II ♦ - Example H considers the disk studies in example I 
at the time that the operating conditions have reached those indicated 
by curve I-I of figure 6. The elastic calculations are made by the 
method of reference 1 modified in accordance with the changes made 
in example I. The essential parts of the computation are shown in 
table II, which is abridged from the complete calculation. Column 25a 
lists the values of [P' n ] that were found as the final values 
of P’ n in example I. Plastic flow occurs at stations 17, 18, 19, 
and b, and calculated true stresses when this plastio flow is 
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considered are listed in table II. The stresses obtained as a result 
of this computation are plotted in figure 9. The elastic stresses 
obtained without considering the plastic flow that occurred 
previously are also plotted for comparison in figure 9. 

Example III. - Example HI continues the cycle analyzed in 
examples I and II, at the conditions of curve IH of figure 6. 

Table IH gives the essential parts of the calculation for this 
example, which is similar in procedure to table I. The value of 
[P 'nj, column 25a, however, is the total of the values of P’ n 
obtained from examples I and H. In this example, plastic flow 
occurs only at stations 17 and 18. The results of this computation, 
together with the elastic-stress curves found without consideration 
of previous plastic flow, are shown in figure 10. 

Example IV. - The steady-state operating conditions represented 
by curve IV of figure 6 are treated in example TV. The essential 
calculations shown in table IY(a) were made similarly to those in 
table III except that [P' n ] in column 25a is now the sum of the 
values of P' n from examples I, H, and IH. Because no values 
of a e,n exceed those of Oj } n , no plastic flow occurs and the 
stress values of table IV(a) are the true stresses at the beginning 
of steady-state operation. However, as parts of the disk are at 
elevated temperature, significant creep can occur at steady load 
at stations 16, 17, and 18, where the stresses are sufficiently 
high. Table IV(b) shows the calculations of creep. Column 63 lists 
the creep rate Cq (in./(in.)(hr)), and column 64 the creep 
increment T n for the 5 -hour running period. Columns 65 and 66 
give the computed values of 5 r ^ n and S-fc,n, respectively. The 
computation then proceeds in a manner similar to bhe plastic-flow 
calculations, as indicated in the column headings. The values for 
stress obtained indicate that, for small values, creep has only a 
slight effect on the stresses. Figure 11 shows the stress distribu- 
tions at the beginning and the end of the steady- state running 
period, together with the elastic stresses obtained without considering 
either creep or previous plastic flow. 

Example 7. - The conditions of example Y represent one of the 
conditions through which the turbine disk is assumed to pass during 
the stopping period. The abridged elastic calculations are given 
in table Y. Values listed now represent the accumulated effect of 
plastic flow [P' n 3 plus the additional effect of the creep repre- 
sented by [Q' n l [Q' n ] is the s® 111 ® as the Q' n computed in 
example IY. AH values of a e ^ri are less than the corresponding 
values of <7y,nJ therefore no plastic flow occurs. The results of 
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the calculation are plotted in figure 12, together with the elastic 
stresses computed without considering previous plastic flow or creep. 

Examgle 71. - Example YI is the computation of the stress dis- 
tribution at the temperature distribution assumed to be present 
shortly after the wheel has stopped turning. The essential partB 
of the calculation are shown in table VI. Because no flow was found 
in example V, the [P' n ] + [Q’ n ] term will be the same in this 
example as it was in example V. Plastic flow occurs at station b. 

The resulting stresses are plotted in figure 13, together with the 
elastic stresses computed without consideration of previous plastic 
flow and creep. 

Example VH. - Example VII is the computation of the stress 
distribution in the dish after the temperature has become uniform 
at the ambient temperature (assumed to be 70° F) throughout the disk. 
These stresses are therefore the residual stresses in the disk 
resulting from the flow occurring during the complete operating 
cycle. The abridged calculations given in table VII indicate that 
plastic flow occurs at station b and the residual stresses are 
plotted in figure 14. 

Discussion of numerical examples. - The foregoing cycle of stress 
calculations is indicative of the means of obtaining a complete 
analysis of the stress behavior of a turbine disk. Although the 
results plotted in the various figures and summarized in figure 15 
do not represent the exact behavior of any particular turbine disk 
because of the lack of data on the material properties and temper- 
ature gradients, they do give a qualitative picture of the behavior 
of a turbine disk with welded blades. The high residual tensile 
stress at the rim of the wheel provides a plausible explanation of 
the rim cracking that has occurred in such wheels. The compressive 
flow at the rim during starting and the tensile flow on stopping 
result in cyclic flow of the rim material with each start and stop 
and possibly induce cracks. When accurate data are available on 
creep, stress-strain curves, the effect of strain-hardening, and 
temperature distribution, quantitative analyses of disk behavior will 
be available as a guide in future turbine design. 


CONCLUSIONS 

A method for studying the operating stresses in gas-turbine 
disks has been presented that includes consideration of the effect 
of plastic flow and creep on the stress distribution. Results of 
calculations indicate that rim cracking in turbine wheels with 
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welded blade attachments may be caused by alternate compressive and 
tensile plastic flow as the wheel is alternately heated and cooled. 
From the results of the numerical examples presented, it may be 
qualitatively concluded that plastic flow alters the elastic-stress 
distribution markedly and that, if the amount of creep is small, the 
effect on the stress distribution is also small. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 5, 1948. 
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APPENDIX - STRESS CALCULATION FOR DISK 
WITH CENTRAL HOLE 

The calculations given in the section of the report entitled 
"NUMERICAL EXAMPLES" dealt with a dish that was solid at the center 
and had temperature gradients such that the plastic flow was 
confined to the region of the rim. Disks of other types spun under 
different conditions may he subject to plastic flow in other regions. 

One example of such a disk is a parallel-sided disk with a 
central hole spun at a uniform temperature. In this disk plastic 
flow first occurs in the region of the central hole. Such a disk 
spun at a speed great enough to cause some flow near the hole is 
calculated here. 

The essential columns of the elastic calculation are given in 
table TUI (a). Flow is indicated at stations a, 2, 3, and 4. How- 
ever, as flow occurs the stresses farther out in the disk may be 
increased. The quantities depending on disk dimensions, together 
with the first approximation, shown in table VHI(b), are found in 
the manner given in the text. However, when the values of Bj, n 
and B-fc^n (columns 57 and 58) are found for the stations at which 
flow occurs, new values of B^n and B^ n must be computed for 
all other point stations also before a new value of at, a 
(column 59) can be found. These computations are also shown in 
table Vin(b) for the first approximation. Additional approxima- 
tions must be made in the same manner until the correct flow 
increments are found. The stresses so calculated are> plotted In 
figure 16. 

Where large numbers of computations involving plastic flow at 
the center of the disk are to be made, it may be desirable to change 
the finite-difference approach to the problem in such a manner that 
the calculations are made from the outside of the disk toward the 
center instead of from the center towards the rim. This procedure 
has certain disadvantages as a general approach to the problem of 
stresses in disks, particularly in that it requires a greater 
number of significant figures to obtain the same accuracy. However, 
for special applications it may present advantages that outweigh 
the disadvantages in more general problems. 
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-1,627.2 
-3,018.3 
-4,744.9 
rl0,746 
-20,313 
-37,271 
-28.938 
-38,664 

1.0000 

.99999 

1.0000 

.99998 

.99996 

.99997 

.99996 

1.1239 

1.3896 

1.4309 

1.7019 

1.7971 

1.0120 

1.7941 

1,4231 

1.4303 

1.6013 

1.6057 

1.0000 

1.0000 

1.0000 

.99999 

.99996 

.99999 

.99999 

1.0591 

1.1643 

1.0916 

1.3346 

1,4640 

1.4941 

1.5149 

1.3080 

1.3497 

1.3843 

1.4011 

tm 

0 

-34.786 
-76.630 
-170.44 
-306.64 
-441.97 
-641.33 
-B.173.8 
-4,104.0 
-T, 367.8 
-9,457.7 
-11,679 
-13,403 
-16,968 
-16,800 
-20,718 
-26,267 
-31,137 

0 

-13.298 

-30.034 

-96,431 

-178.39 

-266.06 

-609.83 

-1,887.3 

-9,470.6 

-4,813.8 

-6,997.5 

-10.469 

-14,344 

-88,916 

-44,088 

-79,644 

-100,070 

-139,900 

-181,850 

-213,300 

i L 

t s . 
a s a 

i w 

3 o 

& * 

* 

u 

D 

U 

87,766 
27,719 
87,578 
87,676 
27,447 
27,298 
86,893 
29,076 
31,600 
36,288 
- 37,779 
50,300 
56,689 
33,043 
23,000 
10,702 
19,412 
13,374 
9,693 

ix£2S 

87,766 

87,740 

87,718 

87,669 

87,369 

27,488 

87,845 

28,100 

39,659 

30,938 

31,493 

39,958 

as, 103 

13,131 

-4,504 

-41,628 

-67,593 

-101,090 

-141,780 

73,500 

73,500 

73,500 

73,600 

73,500 

73.600 
73,500 

78.600 

73.600 
73,600 
73,600 
73,600 
73,600 
73, *00 
73,800 
73,500 

72.000 
70,300 

67.000 

Il.lWl 

87,786 
87,730 
87,898 
87,618 
27,816 
87,391 
87,071 
20,404 
30,636 
33,886 
36,061 
34,883 
38,486 
89,364 
37,639 
33,490 
77,180 
106,400 
149,760 
JMS 
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TABU I. - CALCULATION CP 3THSB3 DI3THIBT7TIOR FOR EXAMPLE I - Oontinned 

(b) Cana tint a determined by dielc goor»trj to 
bo uiM in nil plaatio oalonlatlone 



42 

43 

44 

43 

n 

(9)/<l) 

(3>/(l>n-l 

1 + (48 } 

1 - (43) 

17 

0.016162 

0.016883 

1.0168 


XU 

.014700 

.016138 

1.0147 


19 

.014888 

.014706 

1.0143 


b 

• 013889 

.014886 

1.0139 



(o) Plaetia-atreee onloolntitm 



46 

48a 

47 

48 

49 

60 

51 

62 

63 

Approx- 

imation 

n 

*p,n 

(graph 

reading, 

t lg.7) 

* P.n 

(estimate ) 

*r,n> 

[^ 27 ^ 1 

£8(60i-(61)J 

*t,n» 

(48ni 1 
(6B)z2 

CS(61 )-(6u 33 

(47)x(42) 

(47) n . 1 x(43) 

(48)1(44) 

(48) n _ 1 x(45) 

p v 

(49)+(80) 

-(B1)+(B2) 


17 

20xl0“« 

aoiio' 6 

13.081x10-4 

-'19.057x10*6 

0.19729x10-® 

0.00000x10-4 

-19.966xl0-« 

0.000x10-4 

20.153x10"® 

1 

10 

870 

670 

398.07 

“661.16 

6,0099 

0.197 29 

-676.94 

-19.369 

662.39 


19 

1860 

I860 

1014.8 

-1847.1 

14.489 

6.8099 

-1873.6 

-636.36 

1237.4 


b 

3960 

3080 

2038.6 

-3989,1 

28.864 

14.489 

-4014.1 

-1820.7 

2236.8 


17 

80 

20 

12.362 

-19.769 

0.19034 

0 

-20.069 

0 

20.249 

a 

18 

090 

686 

402.14 

-681.31 

3.9139 

.19034 

-691.33 

-19.460 

677.97 


IS 

IS 30 

1S0O 

1040*5 

-1360.7 

14.66w 

D.axnw 

-1867.3 

-671.B9 

1M5U.U 


b 

4630 

4300 

2312.1 

-4296.8 

32.115 

14.868 

-4565. 6. 

-1834.1 

8668.4 


17 

80 

80 

12.344 

-19.762 

0.19007 

0 

-20.002 

0 

20.868 

3 

18 

370 

870 

393.03 

-068.43 

3.7802 

.19007 

-670.23 

-10.463 

662.74 


19 

1BS0 

I860 

1033.5 

-1866.8 

14.807 

8.7802 

-1882.3 

-666.63 

1340.3 


b 

1900 

4100 

2818.3 

-4096.2 

30.810 

14.807 

-4162.1 

-1889.3 

2368.4 


17 

20 

20 

12.688 

-19,760 

0.19023 

0 

-20.060 

0 

80.260 

4 

1U 

ovu 

bbO 

098.99 

-676.50 

3.B075 

.19023 

-686.30 

-19.461 

072.90 


19 

I860 

1860 

1037.2 

-1866.7 

14,817 

6.8873 , 

-1883.2 

-666.42 

1236,6 


b 

8900 

4000 

2166.0 

-3996.8 

29,946 

14.817 

-4061.3 

-1889.2 

8266.2 


17 

80 

20 

12.568 

-10.769 

0.19028 

0 

-20.059 

0 

20.249 

5 

18 

680 

080 

399.27 

-678.33 

6.8717 

.19028 

-686.27 

-19.460 

672.87 


19 

I860 

1880 

1037.5 

-1866.7 

14,822 

6.8717 

-1883.8 

-660.38 

1236.8 


b 

4200 

4016 

2160.3 

-4011.0 

30,004 

14.828 

-4066.8 

-1829.8 

8882.4 


17 

20 

80 

12.637 

-19.760 

0.19026 

0 

-80.060 

0 

20.260 

0 

18 

060 

680 

399.26 

-676.33 

6.8714 

.19026 

-686.27 

-19.401 

678.87 


19 

I860 

1880 

1037.3 

-1886.7 

14.822 

3 .8714 

-1808.2 

-666.38 

1266.5 


b 

4018 

4018 

2180.8 

-4011.0 

30.011 

14.882 

-4068.8 

-1829.8 

2282.4 




to 

\JT 
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a T hl i rmlvto of i# *1'° substitubii for tha original indue of O" t,a ul0< * in table 1(a) to compute plastlo etresa for stations a to 16. 
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table n. - AEanxjm vaujes from calotiaiiob op strem distrihjtiob for example n 



i 

£ — 

6 

6 

7 


SB* 

isb 

36 

39 

34 

60 

01 

n 

r n 

*n 


a n 

4 *n 

H 'n 

[py + W t ] 

H' n -ip-j -wy 

°Y,n 


H'n-CPtf-CQ y 
-PU 

°"r,n 


A. 

o 

3 

4 

B 

a 

7 

8 
9 

10 
1 1 

1 2 

13 

14 

16 

16 

17 

10 

19 

b 

0.6000 

AQCrv 

]7flOO 

1,0000 

1.2600 

1.5000 
2,0000 

3.0000 
4*0000 

5.0000 

R.fUYYl 

6.0000 
8.6000 

7.0000 

7.5000 

8.0000 
0.2600 
0.6000 
8.7B00 
9,0000 

4.3750 

a 

n«l ww 

4.3750 
4.3780 
4.3760 
4.3750 
4.3760 
3.0400 
3,2760 
2.6000 
C» . ^TWt 

2.2100 
a . 1600 

2.1650 

2.7000 

2.7000 

2,5600 

8.3800 

8.1460 

1.9100 

29*700x10® 

OO 'T/in 

29 ! 700 

89.700 

29.700 

89.700 

80.700 

20.700 

29.700 

29.700 
QQ-Ann 
29 ! 600 

29.400 
89,800 

28.700 
27.900 
27.300 

20.400 

34.700 

20.700 

8. 6060x1 0*« 

o nrjtn 

eisooo 

0.6060 

6.5060 

8.5060 
B, 3050 
0.5060 
0,6070 
0.6100 
fl.MVl 
0.6620 
0,6130 
8.7000 
6.0410 
9.0680 
9.8260 
9.4200 
9.6670 
9.9490 

130 

1 

4.»JV 

130 

130 

130 

130 

130 

130 

131 
138 
147 
163 
197 
851 
339 
401 
379 
700 
046 

1030 

A — 1 

yj juu - 

0 
0 
0 
0 
0 
0 

8,6400 

31.060 

no Ctnr\ 

(u.uiW 

136.36 

204.03 

483.94 

613.40 

1364.6 

900.14 

1252.2 

1693.1 

QfiFiA . % 

0 xIO** 
0 
0 
0 
0 
0 
0 
0 
0 
rt 
V 

0 

0 

0 

0 

0 

80.260 
372.87 
1233,3 
gang . a 

G xiG” 6 

0 
0 
0 
0 
0 
0 

0.6400 

61.060 

io,S70 

158,30 

284,03 

403.94 

813.40 

1334.0 

959,69 

579.33 

530.60 

- 224*10 

33,093 
36,037 
36,973 
35,014 
36,613 
33,370 
34,746 
57,840 
39,977 
43,774 
45 , 415 
46,390 
43,687 
39,276 
83,333 
19,076 
15,833 
13,469 
10,192 

n ^ Q£2 

33,093 

36,062 

33,033 

33,944 

33,821 

36,675 

35,296 

33,016 

37,657 

36,808 

37,549 

34,173 

86,320 

12,677 

-18,109 

-47,307 

-70,310 

-00,607 

-01,666 

=51,549 

0 tHQ = * 

0 
0 
0 
0 
0 
0 

8.6400 

61.030 

7G.S70 

168,38 

004.03 

406.94 
013.40 

1334.6 

036.94 
300,00 
251.64 

=10 5. 13 “ 

35,305 

35,689 

33,333 

36,403 

36,806 

34,962 

34,330 

33,701 

39,451 

43,149 

45,720 

45,567 

42,887 

36,643 

23,934 

19,296 

16,301 

13,031 

9,908 

rr nflo 

f t \rum. 

36,303 
33,661 
35,687 
33,633 
33,413 
36,237 
34,007 
36,783 
37,003 
37,700 
ov,yu3 
33,801 
33,710 
12,060 
-12,673 
-47,941 
-67,610 
-73,003 
-72 , ioe 
-66,196 


TABCE XU. - AHUDOBD TALDE3 FROM OALCtfLATIOE CP STRESS DISTRIBUTION Fffit EXAMPLE III 



i 

s_ 

8 

3 

7 

25 

95ft 

26b 

3© 

39 

53 

55 

31 


r n 

**n 

“n 

a n 

«n 

H'n 

CP’j+Ci'j 

H' -C*U -ran 

n " “ -■* n 

O" 

A’,U 


tT 1 _f m 1 r~ 

“ D 

-P’n 

~r,n 


A 

2 

3 

4 

3 

6 

7 

a 

9 

10 

11 

IS 

13 

14 

15 

13 

i n 
18 
ie 

b 

0.5000 

.6260 

.7500 

1.0000 

1,9500 

1.5000 
2.0000 

3.0000 

4.0000 

5.0000 

5.3000 
0.0000 

6.3000 

7.0000 

7.5000 

e.oooo 

n nmris \ 

W(A LW 

8.5000 
0.7500 

9.0000 

4.3750 

4.5750 

4.3780 

4.3750 

4.3750 

4.3750 
4.3750 
5.0400 
3.8750 
3,3800 
2.3720 
2.2100 
2.1500 
2.1660 

2.7000 

8.7000 

O ■ 

2 ! 3600 
2.1450 
1.9100 

20.9001X1® 

28,000 

88.800 

20.000 

88.600 

20.000 

20.000 

28.800 

20.700 

88.700 
2S.600 

08.400 
20.200 
27,900 

27.300 

26.400 

O M 

iV.SOO 

20.400 

19.300 

0.0860x10”® 

8.8230 

0.0260 

0.6230 

0.0250 

0.8960 

8.8260 

8 .0 £oO 

8.8330 

0.0620 

8.0090 

6.9310 

0.9950 

9.0930 

0.2290 

9.4900 

0 *AriA 

b!6790 

9.8410 

10.039 

330 

330 

330 

330 

330 

330 

330 

551 

336 

359 

339 

.395 

436 

493 

5B1 

700 

rnr» 

1 1 V 

008 

963 

1000 

0 xltf® 

0 
0 
0 
0 
0 

9.DQO0 

48.000 

160.30 

100.60 

947.70 

394.60 

307.00 

631.90 

1232.0 
799.50 
949.80 

1133.6 

1334.1 

0 xicf 8 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

141.19 

072,12 

1345.5 

9183.4 

0 ya<r® 

0 
0 
0 
0 
0 

0.5000 

43.000 

160.50 

150.00 
247.70 
394.60 

507.00 
031,90 

1232.0 

630.51 
77.60 

-209,90 

-809.30 

44,433 

44,361 

44,277 

44,066 

43,002 

43,480 

42.655 

48,497 

48,599 

52,198 

34,770 

54,069 

50,060 

44,243 

29,200 

21,358 

10,041 

14,659 

11,007 

9,339 

44,435 

44,403 

44,339 

44,230 

44,079 

43,884 

43.301 

44,159 

44,051 

43,401 

41,105 

35,838 

25,236 

9,000 

-15,119 

-44,820 

-58,300 

-66,752 

-44,468 

-21,474 

0 xio" 6 

0 
0 
0 
0 
0 

9.5000 

46.800 

130.50 

160.60 

347.70 

394.60 

587,80 

851.90 

1238,0 

633.34 
91.355 

200.34 
-009,30 

44,410 
44,344 
44,260 
44,049 
43,786 
43,463 
42.336 
45,' 478 
48,377 
59,172 
54,749 
64,032 
50,036 
44,812 
09,178 
21,320 
18,028 
14,666 
11,007 
9,330 

44,418 

44,306 

44,342 

44,221 

44,063 

43,067 

45,364 

44)141 

44,642 

43,459 

41,102 

35,533 

96,218 

9,834 

-15,136 

-44,031 

-57,901 

-38,325 

-44,470 

-21,403 


hO 
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TJtBLE 17. 


OALCULATIOH OP 8»S3fl DI3MIE1P7105 PCE SXAKFUl XV 


K5 

OD 


(a) Abridged values 



1 

8 

6 

« 

7 

85 

26a 

86b 

38 

39 

72 

70 

70 

n 

*n 

V 


a n 


H* 

n 



""r.n 

•*,n 

“’n-U'iJ-MU 

-Q'n 


°*,n 

a 

0.6000 

4.3760 

37.600x10® 

9.1470x10“® 

630 

0 xlO" 6 

o xicr® 


48,000 

46,006 


46,006 

46,005 

a 

.6860 

4.3780 

87.600 

9.1470 

630 

0 

0 

0 xlO" 6 

46,938 

46,974 

0 xlO" 6 

46,931 

48,973 

8 

.7600 

4.3760 

87.600 

9.1470 

630 

0 

0 

0 

46,848 

46,930 

0 

46,847 

45,920 

4 

1.0000 

4.3700 

87.800 

9.1470 

630 

0 

0 

0 

46,637 

46,809 

0 

46,636 

46,800 

6 

1.8600 

4.3780 

87.800 

9.1470 

630 

0 

0 

0 

46,373 

46,660 

0 

48,378 

46,649 

6 

1.6000 

4.3?60 

67.800 

9.1470 

630 

0 

a 

0 

48,061 

46,466 

0 

45,060 

46,464 

7 

8.0000 

4.3760 

87.600 

9.1490 

631 

10.200 

0 

10.800 

44,194 

44,698 

10.300 

44,193 

44,697 

0 

3.0000 

3.8400 

87.800 

9.1680 

537 

59.700 

0 

69.700 

46,969 

44,419 

69.700 

46,968 

44,418 

9 

4.0000 

3.8780 

87.600 

9.1840 

853 

161.00 

0 

161.00 

49,449 

48,607 

161.00 

49,448 

48,606 

10 

6.0000 

8.6800 

87 .'300 

9.8380 

587 

343.90 

0 

343,90 

62,329 

37,156 

343,90 

68,328 

37,184 

11 

6.6000 

8.3780 

87.100 

9.8880 

614 

276,40 

0 

876.40 

64,197 

32,627 

276.40 

64,198 

32,626 

IB 

6.0000 

0.8100 

86.800 

9.3380 

649 

361.30 

e 

361.30 

62,668 

26,030 

361.30 

62,666 

83,026 

13 

6.6000 

8.1600 

86.800 

9.4000 

693 

469.30 

0 

469.30 

47,980 

14,400 

459.30 

47,918 

14,399 

14 

7.0000 

8.1660 

88.900 

9.6000 

730 

608.30 

0 

606.30 

41,617 

927 

606.30 

41,616 

'227 

18 

7.6000 

8.7000 

86.100 

9.6100 

819 

748.60 

0 

746.60 

86,780 

-19,803 

746.60 

86,778 

-19,204 

16 

B.0000 

8.7000 

83. BOO 

9.7480 

906 

961.30 

0 

961.30 

19,146 

-36,637 

961.30 

19,146 

-39,628 

17 

8.8600 

8.5600 

28.700 

9,9870 

954 

706.70 

166.17 

639.63 

10,107 

-46,488 

638.64 

16,108 

•45,443 

ie 

8.6000 

8.3800 

81.800 

9.9120 

1007 

463.80 

868.44 

-404.64 

13,200 

-SB, 446 

-406.17 

13,199 

-flBiilfl 

19 

8.7600 

2.1480 

19.800 

L0.007 

1006 

686.60 

1334.0 

-648.40 

11,008 

-16,684 

-647.94 

11,003 

-16 *,333 

b 

9.0000 

1.9100 

17.800 

10.109 

1130 

766.00 

8163.4 

-1407.4 

9,340 

10,470 

-1406.6 

9,334 

10-470 


(b) 0*1 oulet Ion of offset of oroop on final rtraaeea 



63 

64 

66 

66 

67 

M 

69 

70 

71 

n 

°n 

r tt4 

(63 )x6 

^r,n» 

■jJn^[8(3eM39H 

^t,n» 

7 rf^r- x[b(m )-(»}] 

(66)1(42) 

(66) n _ 1 x(4B) 

(66)1(44) 

{66) n . 1 x(46) 

(07)+(68 )- 
(69 iei7o i 

17 

0.2x10-® 

1.0x10-® 

0.70226x10** 

-0.96766x10“® 

0.010641x10"® 

0,00000x10"® 

-0. 98237x1 (T® 

-O.OOOOOxlO" 8 

0.99301X10" 6 

10 

• 3 

1*5 

1.0848 

-1.4596 

.016963 

.010641 

-1.4608 

-.96300 

.63439 

19 

*2 

1.0 

.#024' 

-.91800 

.011464 

.016963 

-.93113 

-1.4184 

-.43986 

b 

0 

0 

0 

0 

0 

.011464 

0 

-.90488 




75 

73 

74 

76 

76 

77 

78 

79 

B 

H'-tP'J-CQ'J 

"0'n» 

(26b ) -(71 ) 

“n* 

t(72)xfl0) + 
il4)x(81IH (86] 

„ H ’n* 

I78)x<8) + 

(14 )x(20 )] ♦ (26) 

8 r a , 

(27 )x(75) n _ 1 + 
(28)x(76^ n -l 

B t,n* 

(29)1(75 Jn-i* 

or 

t,a» 

r^, b -( 7B >b3 

+(M) b 

'’r.n* 

(33)x(77)f(7B) 

(34)1(775+ (76) 

17 

638.54x10"® 

-2469.1 

-18,076 

-62,373 

-103,190 

46,003 

16,105 

-46,443 

18 

-406.17 

-8847.2 

7,607.6 

-59,834 

-91,838 


13,199 

-32,418 

19 

-647.94 

-2306.8 

11,761 

-69,620 

-76,973 


11,003 

-16,636 

b 

-1406.5 

-8830.7 

83,906 

-60,379 

-60,965 


9,334 

10,470 




0£6 
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! TABLE V, - ABRIDGED VALUES TECS OALCOLATIOB OP STRESS DI3TOUUTI0H FCH EXAMPLE V 



i 

8 

s 

6 

7 

26 

85a 

25b 

38 

39 

54 


81 

Q 

r n 

\ 

E a 

a n 


H ’n 


BBSS! 

El 

m 


l 

°t.» 

a 

& 

0,5000 

.6930 

4.5750 

4.5780 

88.800x10® 

88.800 

8 . 8280x10"® 
8.8500 

331 

358 


0 Xlo -6 

9.49 xlO-® 

83,288 

28,263 

£2,205 

22,023 

9.49 xlO*® 

28,883 

22,863 

£2,205 

22,023 

5 

.7500 

4.5750 

38.800 

8.8510 

333 

9.16 

0 

9.16 

EE, 181 

21,811 

9.16 

82,181 

31^811 

4 

1*0000 

4.5750 

28.800 

8.8350 

334 

9.60 

0 

9.60 

£2,027 

81,645 

9. BO 

88,087 

21^646 

6 

1.8500 

4,3730 

28,700 

8.8390 

336 

37,30 

0 

37,36 

81,801 

£0,667 

37.36 

81,801 

£0,667 

6 

1.60 CO 

4*5750 

28.700 

8.8440 

341 

88.32 

0 

28.82 

21,501 

80,005 

£6,23 

81,601 

20 j 005 

7 

2,0000 

4.5760 

28.700 

B.8680 

550 

84.50 

0 

84.50 

20,744 

18,160 

84 .50 

20,744 

18 j 160 

0 

5.0000 

3. 8400 

28.500 

8.8970 

374 

288.70 

0 

826.70 

81,090 

14,149 

226,70 

81,090 

14', 149 

9 

4.0000 

3,9760 

88.500 

8.9550 

409 

354.76 

0 

554.78 

80,904 

8,300 

334.78 

80,904 

8^500 

10 

6.0000 

2.6800 

98.100 

0.0840 

453 

486.09 

0 

486.09 

80,680 

1,010 

486.09 

80,680 

1^010 

11 

6.5000 

8.5720 

88.000 

9.0650 

479 

£34,27 

0 

£34,27 

20,586 

-3,348 

834.87 

80,688 

-3^342 

12 

0.0000 

8.2100 

87,800 

e.nso 

608 

£87.75 

0 

867.76 

19,144 

-8,784 

287.76 

19,144 

-6,724 

13 

5.5000 

2.1600 

87.600 

9.1620 

639 

309.48 

0 

309.48 

18,689 

-14,887 

309.42 

16,589 

-14^087 

14 

7.0000 

2.1550 

87.400 

6.2140 

572 

538.09 

0 

339.09 

13,400 

-81,658 

538.09 

13,400 

“21^668 

IS 

7,5000 

8.7000 

87.100 

9,2780 

008 

365.97 

0 

366.97 

7,748 

-89,867 

306.87 

7,743 

-29^807 

16 

8.0000 

2,7000 

86.600 

9.5350 

645 

391.74 

0 

391.74 

4,366 

-37,781 

391.74 

4,366 

-37^781 

17 

8.B500 

8.5600 

36,700 

9.3630 

666 

807.07 

167.10 

40.840 

8,879 

-37,439 

40,840 

8,879 

-37^429 

18 

8.6000 

2,5800 

26.600 

9.3990 

687 

880.03 

858.97 

-638.97 

1,676 

-19,393 

-638.97 

1,675 

-19^393 

19 

8.7600 

8.1450 

26.500 

9.4490 

708 

339.78 

1535.6 

-1100.5 

1,171 

9,033 

-1100.6 

1,171 

0,633 

b 

9.0000 

1.9100 

86.100 

9.4990 

730 

244.58 

8108,6 

-1910.5 

1,764 

58,858 

-1918.8 

1,764 

68i868 


TABU? VI. - ABRIDGED VALUES FHOt OALOOLATIOE OP STRESS DIBTHIBDTIOI FOR EXAMPLE VI 



1 

8 

- -B 

6 

7 

86 

26a 

§£Tb 

36 

39 

64 

55 — 

61 

n 

P 

a 

h 

n 

5 n 

a 

n 


H, b 

CP'JKQ'J 

H’ -CP’J-W 1 

n n n 

^r.n 

^,11 

HV-cpy- [tl y 
- P'a 

cr 

r,n 

0“ 

t,n 

a 

2 

3 

4 

5 

6 
7 
S 
9 

10 

11 

18 

13 

14 

15 

16 

17 

18 
19 

_b 

0.6000 

.6250 

.7500 

1.0000 

1.3600 

1.5000 
2.0000 

3.0000 

4.0000 

5.0000 

5.5000 

6.0000 

6.5000 

7.0000 
7,6000 

8.0000 
8.8500 

6.5000 
6.7500 
9.0000 

4.3750 

4.3760 

4.3760 

4.3760 

4.3760 

4.5760 

4.3760 

3.8400 

3.8760 

8.6800 

8.3780 

8.8100 

8.1600 

2.1560 

2.7000 

8.7000 
2.5600 
3.5800 
8.1460 
1.9100 

29.300x10® 

89.300 

89.300 

89.300 

29.200 

89.800 
89.600 
89.100 
89.100 
89.000 

89.000 

28.900 

88.900 
88.900 

88.000 

88.800 
88. BOO 
88. BOO 
88.800 
88.800 

8.6760x10-® 

8.6770 

8,6790 

8.6840 

8.6880 

8.6930 

8.7010 

8.7100 

8.7370 

8.7660 

B.7640 

8.7730 

8.7820 

0.7910 

8.7990 

8,8090 

8.8140 

£.0170 

8.8880 

8.8260 

836 

837 

838 
241 
244 
247 
862 
263 
874 
886 
891 
297 
308 
308 
313 
319 
388 
384 
387 
330 

B. 91x10-® 
9.16 
87.20 

87.03 
27.30 
43.48 

100.44 

100.84 

110.88 

46.11 

66.88 

46.68 

66.46 

46.46 
66. 9B 

88.04 
18.60 
88.09 
27.79 

o Xlo* 6 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

167.10 

858.97 

1333.5 

8168.6 

8.91x10-® 

9.15 

27.80 

27.03 

27,30 

43.48 

100.44 

100.84 

110.38 

46.11 

53.90 

46.68 

85.46 

46.46 
86.96 

-139.18 

-840.37 

-1306.4 

-8134.7 

3,616 

3,698 

3,537 

3,364 

3,140 

2,697 

8,417 

1,785 

933 

2 

-643 

-1,108 

-1,639 

-8,140 

-3,163 

-2,676 

-3,033 

-3,000 

-2,197 

0 

3,616 

3,376 

3,100 

8,629 

1,960 

1,403 

657 

-1,301 

-3,144 

-5,260 

-6,093 

-7,186 

-8,043 

-9,160 

-9,819 

-10.8B9 

-6,736 

17,132 

63,634 

118,600 

8.91x10-® 

9.13 

27.20 

37.03 

37.30 

46. 45 
100.44 
100. 34 
110.88 

46.11 

36.26 

46.68 

66.46 
46.43 
55.98 

-139.13 

-840.37 

-1306.4 

-1461.0 

3,760 

3,726 

3,678 

3,498 

3,274 

3,030 

8,351 

1,876 

1,106 

307 

-317 

-868 

-1,396 

-1,900 

-1,974 

-8,488 

-2,884 

-2,7B6 

-1,969 

0 

3,750 

3,609 

3,894 

8,668 

8,090 

1,536 

690 

-1,160 

-2,990 

-0,087 

-6,907 

-6,991 

-7,843 

-8,947 

-9,031 

-10,708 

-6,664 

17,329 

63,740 

93,678 
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TAHIZ VU. - A EE IDG ED VALDES FROM OALCCLATIOH OP STRESS DISTRIBOTIOS POE EXAMPLE VII 


OJ 

o 



i 

2 

6 

6 

7 

26 

36* 

26b 

38 

39 

64 

60 

01 

n 

r n 


% 

a n 

6T n 

H’n 

trynug 


^r.n 

«t, Q 

H'n-CP'J-^y 

- P'n 

“"r.n 

,Q 

I 

1 

0.0000 
• 6250 
.7000 
1.0000 

1.2500 
1.5000 
2,0000 

5.0000 
■4.0000 

6.0000 
6.6000 
6.0000 
6.6000 

7.0000 

7.6000 

8.0000 

8.2500 

8.6000 
0.7500 
9.0000 

4.5780 

4.3750 

4.3750 

4.3750 

4.3750 

4.3750 

4.3760 

3.0400 

3.8760 

8,6800 

8.3780 

8.3100 

3.1800 

8.1560 

8.7000 

8.7000 

2.5600 

2.5600 

2.1450 

1.9100 

to 

0 

• 

8 

•P 

4 

1 

f 

o 

3 

CO 

*P 

■P 

s 

O 

o 

u 

1' 

a 

o 

o 

o 

■p 

! 

u 

0 xio -8 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

167.16 

858.97 

1333.6 

1488.8 

0 xlCf* 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

-107.16 

-858.97 

-1333.6 

-1488.8 

-2,408 

-2,468 

-8,468 

-8,468 

-8,468 

-8,408 

-8,408 

-8,779 

-3,183 

-3,779 

-4,802 

-4,437 

-4,476 

-4,432 

-3,51fl 

-3,615 

-3,044 

-3,378 

-8,292 

0 

-8,468 

-8,468 

-8,408 

-3,460 

-8,408 

-3,468 

-3,468 

-8,014 

-8,849 

-3,193 

-3,438 

-3,604 

-3,706 

-3,768 

-3,470 

-3,479 

1,368 

26,518 

68,023 

107,180 

0 xio* 6 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

-167.10 

-B58.97 

-1333.5 

-947.81 

-2,556 

-3,560 

-2,356 

-2,356 

-8,356 

-2,366 

-8,366 

-2,653 

-3,039 

-3,608 

-4,018 

-4,237 

-4,873 

-4,833 

-3,369 

-3,366 

-3,476 

-3,199 

-2,093 

0 

-2,360 

-3,366 

-8,366 

-9,366 

-8,366 

-8,366 

-8,366 

-8,496 

-2,720 

-3,049 

-3,873 

-3,448 

-3,337 

-3,692 

-3,318 

-3,331 

1,519 

80,677 

66,196 

91,169 


TABLE TUI. - QALOULATIOH OP STRESS DISISIBOTIOS PCH PARALLEL-SIDED DISK 
(a) Abridged taIuob 



ni 

8 

5 

5 

7 

26 

86a 

25b 

30 

39 


41 

.. 40. 

61 

n 

n 

h n 

E n 

a 

n 

AT„ 

H* 

n 



°"r,n 

IBS 


°"e,n 

ES 

mm 

a 

0.5000 









176,040 





£ 

. 6 £50 








0 


176,040 

0 

134,620 

3 

.7500 








31,050 

142,320 


129,620 

24,952 

127,160 

4 

1.0000 
1 ,5000 



b 





48,400 

124,860 


109,030 

41,519 

124,510 

6 


d 

g 





04,873 

106,350 


92,002 

50,031 

111,530 

6 

2.0000 
2. 5000 

o 

o 

8 

£ 

o 

o 

o 

o 

74,394 

91,900 


84,518 

72,549 

93,977 

7 

§ 

§ 

-p 



43 

75,607 

86,673 

8 

61,704 

74,625 

B7.813 

8 

3.0000 

3.5000 

o 

• 

M 

* 

4 

4 

4 

74,187 

B3.266 


79,119 

73,5a© 

B4.006 

© 

rH 


CD 

-y 

£ 


43 

71,369 

80,278 

8 

H 

76,215 

70,976 

80,803 

10 

-P 

-P 

■P 

3 

a 

a 

3 

57,553 

77,268 

72,896 

67,282 

77,669 

11 

4.5000 

6.0000 

fi.fioon 

•3 


d 

■P 

a 

4* 

a 

45 

a 

4* 

68,892 

74,060 

•P 

09,156 

62,700 

74,360 

1ft 

■P 

•P 

x> 

a 

57,471 

70,579 

d 

66,024 

87,332 

70,845 

15 

3 

a 

3 

s 

8 

o 

o 

51,322 

66,782 

■P 

50,551 

51,223 

67,000 

1 4 

fi.oom 


■p 

.3 

a 





44,473 

02,660 

3 

55.8E7 

44,401 ■ 

62,647 

15 


C 

ri 





36,934 

68,163 

■2 

a 

50,979 

36,883 

68,339 

1ft 

7.0000 
7.5000 

8.0000 

8 

§ 






38,711 

63,317 

46,219 

28,078 

63,476 

17 








19,813 

48,103 

8 

41,873 

19,793 

48,249 

b 








10,248 

0 

42,018 

36,662 

38,433 

30,562 

10,333 

0 

43,664 

36,919 
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TABLE VXXI* - OALCDLATIOH OP STRESS DI3TRIHJTI0H FOR PARALL3L -3 U>ED DISK - Concluded 
{b) Caloolation fop first approximation of plaitia-atrosa distribution 


(9)/(l) ( 9 )/ ( 1 ) 


l+(42) 


l-(43) 


0,10000 

.06333 




o, 12500 
.10000 


1.1000 

1.0833 


0,87600 

.90000 


tt p,n P.n 

(graph (oatimats) 
roftding, 
fig. 7) 


1800x10-*] 1800x10 
6 BO | 660 

300 
0 


2(38)-(39)l 


-soo.ooxio-e 

-168.06 

-38.69 

0 


(4i)xe 

[ 2{30)-(30)I 


1800.00x10-6 

657.66 

276.95 

0 


(4?)x(42) 


-16.808x10“® 

-3.2138 

0 


(48)x(44) | (iSJ^x^S 


-112.60x10-6 

-18.806 

-6.4330 


6Bl.S£xl0- 6 

300.02 

0 


) (49)+(60)- 

<61>+<32> 


1576.0x10-6 

483.80 

230.78 


884.37x10 

163.75 

824.34 


H^-cpg-CQyl 


(25b }- (53 ) 


C(54)x(10)+ [<M)x(6) + (37)1(57) 

(14,x(8ia (14)X(20) 3 J 

+ (26 ) +(36) " 1 


-854.37x10-6 
-133. 7B 
-224.34 
0 


2,144. 

117. 

101 . 

-1,878. 

-8,641. 

-3,423. 

-4,218. 

-6,006. 

-5,000. 

-6,697. 

-7,595. 

-8,194. 

-8,994. 

-9,794. 

-10,606 

-11,596 

-18,197 


23,499 

4,476.8 

6,829.7 

-706.88 

-1,060.8 

-1,508.6 

-1,651.8 

-1,799.2 

-8,045.1 

-2,890.8 

-8,534.7 

-8,778.8 

-3,081.9 

-3,206.0 

-3,607.8 

-3,760.1 

-3,992.5 


2.144.1 

5.639.1 
9,836.0 

12,435 

11,484 

9,013,8 

5,667.6 

1,555.5 

-3,588.3 

-8,089.8 

-15,136 

-22,060 

-29,667 

-37,933 

-46,866 

-56,464 

-66,717 


"t.a* *,n* 

(S9)3C(5V ) n _ 1 [<r fb _(57) b ] (33)x(59) 

+(30)x(38) a-1 * ( J 3 ) b +(67) 


£3 
24 

26 
19,838 
17,021 
14,665 
18,309 
9,619.1 
6,619.0 
3,176.3 
-618.0 
-4,670.8 
-9,009.8 
-13,814 
-18,994 
-84,630 
-30,463 


0 

27,816 

43,981 

61,477 


(34)x(69) 


136,510 
135,700 
128, 9B0 
109,960 


138,610 

124,390 

107,930 

96,450 


h«:a tn no 
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Figure 1. - Uniaxial stress-strain curve showing relation 

between stress cr _ and plastic strain € . 

e ,n p,n 
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Strain 

Figure 2. - Uniaxial stress-strain curve illustrating procedure 
used to find correct value of plastic strain. 




Stress 
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Figure 3. - Uniaxial stress-strain curves showing components of 
strain when plastic flow occurs a second time* 
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Figure 4. - Typical deformation-time curve from a constant- 
temperature, constant-load creep test. 




8 


o 


<5 1 2 5 4 5 6 7 

Radial distance, In. 

Figure 6. - Cross section of disk used for numerical examples showing location of point stations. 
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Temperature 


Curve 

Turbine 

Blade loading. 


speed 

°r,b 


(rpm) 

(lb/sq In.) 

I 

8,000 

4596 

II 

10,000 

7062 

III, IV 

11,600 

9339 

If 

V 

e 

1 rri2C 
-L. 1 uu 

vi.vn 

0 

0 


TD 11. 

atiuiLU) f xu« 


gure 6. - Assumed temperature -distribution curves and corresponding turbine speeds 
TCurves I, II, and III are consecutive starting conditions; curve IV represents 
steady-state operation; and curves V, VI, and VII are oonseoutive stopping 
conditions . ) 
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Figure 8. - Stresses in 'turbine disk under conditions of curve I 

of figure 6. 
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Figure 9. - Stresses In turbine disk under conditions of 
curve II of figure 6. 
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Figure 10. - Stresses in turbine disk under conditions of 

curve III of figure 6. 
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Figure 13. - Stresses in turbine disk under conditions of 
curve VI of figure 6. 



Figure 14. - Residual stresses In turbine disk upon completion 

of one operating cycle. 
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Figure 15. - Plastic tangential stresses in turbine disk 
during one running cycle. 
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Figure 16. - Stresses in parallel-sided disk with central hole. 



